Parametric excitation of a dc shuttle current via spontaneous symmetry breaking 
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Abstract 

We investigate theoretically the dynamics of a spatially symmetric shuttle- system subjected to an ac gate voltage. We 
demonstrate that in such a system parametric excitation gives rise to mechanical vibrations when the frequency of the ac signal 
is close to the eigenfrequency of the mechanical subsystem. These mechanical oscillations result in a dc shuttle current in a 
certain direction due to spontaneous symmetry breaking. The direction of the current is defined by the phase shift between 
the ac gate voltage and the parametrically excited mechanical oscillations. Dependance of the shuttle current on the dc gate 
voltage is also analyzed. 



PACS numbers: 73.63.-b,73.23.Hk,78.67.Hc 

Fifteen years ago, a novel form of electron transport 
based on the mechanical vibrations of a metallic nanopar- 
ticle coupled to two electrodes via elastic molecular links 
was proposed by Gorelik et al [1]. Since then, the shuttle 
mechanism of transport has been a subject of intensive 
experimental and theoretical research [2-9] . Different ar- 
chitectures ranging from the original double-junction sys- 
tem to transistor-like configurations comprising a gate 
electrode have been studied [10]. The main feature of the 
orthodox shuttle structures is that a constant potential 
difference, applied between two fixed electrodes, leads to 
a dynamical instability that causes the metal nanoparti- 
cle oscillate. In this regime, the resonator attains a stable 
orbit and a sustained current induced by the voltage drop 
between the electrodes. This current is proportional to 
the frequency oscillation of the nano-oscillator [1]. 

In the present work, we investigate the possibility 
to generate shuttle transport between two electrodes at 
the same electro-chemical potential. In this scheme, we 
demonstrate that, despite the lack of a bias voltage, a 
shuttle dc current can still be detected. This charge 
transport is achieved by applying an ac voltage to a gate 
electrode which controls the electronic population of a 
metallic island and, in this form, also the stiffness of 
the resonator resulting in a parametric mechanical in- 
stability at the resonant frequency. This constitutes a 
new archetype of electron shuttle in which the symmetry 
breaking effect (direction of the shuttle transportation) 
does not rely on the presence of any bias voltage. In the 
phenomena under consideration, the shuttle current is 
controlled by the phase shift between the mechanical vi- 
brations and gate voltage oscillations. We will show that 
in this scenario, two different values of the phase shift, 
which differ from each other by 7r, can correspond to a 
regime of sustained oscillations. The occurrence of these 
values for the phase shift depends, in particular, on the 
initial conditions and, as a result, spontaneous symmetry 
breaking takes place. 



x(t) 




V 6 =VG t +V£ c cofi(a) G t) 



FIG. 1: (Color online) Schematic diagram of the three- 
terminals shuttle device investigated in this work. A quantum 
dot D can oscillate between two metallic leads L, R fixed at 
zero voltage. The dot is capacitevely coupled to a gate elec- 
trode G, to which a voltage V g is applied. Electron tunneling 
takes place between the dot and the leads. Il and Ir are the 
currents between the left and right leads L, R and the dot D. 



In the context of our work, it is worth mentioning that 
parametric excitation in nanoelectromechanical systems 
(NEMS) has been also considered in Refs. [11-14] 

To describe the new shuttling mechanism, we consider 
a system schematically depicted in Fig. 1 where a quan- 
tum dot (D) is connected via elastic links to the left (L) 
and right (R) electrodes. The characteristic dimension 
of the system is d and the mechanical degree of freedom 
of the dot is x(t). The voltage over the left and right 
electrodes is kept null, i.e., Vl = Vr = 0, while a signal 
Vg = Vq + Vq c cosset) is applied to the gate (G). 

In our description, we focus on the single-electron 
regime, where the maximum occupancy of the dot is one 
electron due to Coulomb blockade of tunneling. In addi- 
tion, we also consider the case in which the mechanical 
vibration frequency of the dot is very low in comparison 



to the tunneling rates between the moving quantum dot 
and electrodes. In this framework, the electronic state 
of the dot can be described by the average occupation of 
the dot n(t) and the time evolution of this quantity can 
be written as 



■ =[T L (x)+T R (x)](f-n). 



(1) 



Here, the position-dependent tunneling rates between 
the left/right electrodes and the dot are Tl,r(x) = 
p o6 T^/A w ith A the tunneling length. In addition, / = 
[l + e E D/ k sT^ - 1S Fermi-Dirac thermal distribution, 
where Ejj = e 2 /2Ct(x) — e/3Vc is the electrostatic energy 
of the charged dot, Ct(x) is the total capacitance of the 
system, and f3 « 1 is the transmission coefficient [15]. 
Considering the symmetric case Ct(x) = Ct(—x) and 
small displacements of the dot, x <C <i, we can rewrite 
E D = E - ee 2 x 2 - eV G , with E = e 2 /2C T (0) and 
e & d~ 3 is a positive parameter. 

To consider the mechanical part of the system, we 
model the quantum dot as a single degree-of-freedom os- 
cillator with eigenfrequency uoq subjected to an electro- 
static force Fd induced by the image charges on the leads. 
Having in mind the quasi-adiabatic condition introduced 
before, uoo <C Tl,r, we can make use of the mean field 
approximation Fd = (Fd) = —(dED/dx) = 2ee 2 n(t)x. 
Therefore, the dynamics of the system is described by the 
following dimensionless and coupled ordinary differential 
equations (ODEs), 
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Here, r = uot, = x(t)/d, e = 2e 2 e/raco>Q, a = 

ee 2 d 2 /k B T, v st = (-E + eVg)/k B T, v ac = eV£ c /k B T, 
Co = ujg/ojo, and Q is the oscillator quality factor. 

The average dimensionless current through the system 
can be calculated as: 
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Here, in writing Eq. (3) we use Eq. (2b) and f = 2tt/lo. 
From this expression one can see that a dc current be- 
tween the leads is defined by the correlations between the 



velocity and population of the dot. To find these correla- 
tions one should analyze the dynamical system described 
by Eqs. (2). 

A formal solution of Eq. (2b) is given by, 



n(r) 
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(4) 

Then, by substituting Eq. (4) into Eq. (2a) one can, 
in principle, obtain the exact dynamics of the sys- 
tem by solving the resulting integro- differential equa- 
tion. However, exploiting the smallness of the param- 
eters co>o/ro, e, Q _1 <C 1, a perturbative analysis can be 
employed. In doing so, in the leading order of the param- 
eter cjo/ro, Eq. (2b) reduces to n(r) = /(£, t) and the 
oscillator equation of motion, Eq. (2a), takes the form: 



(5) 



Considering the second harmonic in the Fourier expan- 
sion of /(£, t) in Eq. (5), one can find that the system 
experiences a parametric mechanical instability at the 
driving frequency Co ~ 1. As this takes place, in the 
stationary regime, the time evolution of the dot position 
in the leading order of small parameters e, Q _1 <C 1, 
is given by the expression: 



f(r) = Acos(Cor + x) - 



(6) 



Here, the quantities A = \J2E and x are the stationary 
dimensionless amplitude of the dot and the phase shift 
between the mechanical and gate voltage oscillations, re- 
spectively. The oscillation amplitude, generated by the 
parametric instability, is bounded by nonlinearities and 
it is controlled by the parameter a. 

Turning back to the charge transport, on the accu- 
racy of u;o/ro, the average population of the dot becomes 
n(r) = f(Acos(CoT + x)? r ) using the ansatz in Eq. (6). 
Thus, the expression for the average current of the system 
in Eq. (3) can be recasted as: 



i= d_ f A _JZ sinh(fr|^)|sin( X )) 

2A J_ A CO sh 2 (0 [cosh(a(0) + cosh(6|£(£)| sin( X ))] ' 

(7a) 

with 

a(0 = a £ 2 + "at + bu>£ cos(x) , b = v ac /CbA. (7b) 

Here, |£(£)| is the modulus of the dot velocity as a func- 
tion of its position given by the relation Co 2 \^(t)\ 2 + 
\i(r)\ 2 = 2E. From Eq. (7), for small amplitudes and 
Co ^ 1, the current of the system is given by 



/ = uoQCK^d/ \)A 
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FIG. 2: (Color online) Contour plot for the dimensionless cur- 
rent |/| as a function of the amplitude A and phase x- The 
plot is symmetric with respect to the phase and the sign of the 
current depends on the sign of the phase. The current van- 
ishes at points (black dash lines) A = and/or \ — ±W7r, n — 
0, 1, 2, . . ., current is maximum around x — ±7r/2. 



with 



sinh(^ ac sin(x)) 



cosh(u 3t + y ac cos(x)) + cosh(^ ac sinh(x)) 



(8b) 



From Eqs. (7) and (8), we can conclude that current 
is only attained at a non-zero dot oscillation amplitude 
A while its direction (symmetry breaking) is controlled 
by the phase difference X- This behavior is visulized in 
Fig. 2, which is a contour plot of the dimensionless av- 
erage current as a function of the amplitude and phase. 
From the plot it becomes clear that the sign of the cur- 
rent follows the sign of the phase. We also notice that 
the electronic transport ceases when the amplitude is 
exactly zero A = and/or the phase x = ±^tt, with 
n = 0,1,2..., points highlighted in the phase space in 
dashed black lines. The plot also indicates that the cur- 
rent is maximum around x = ±7r/2. 

Once we have studied the shuttle transport current of 
the system and investigated how it arises from a para- 
metric excitation. We proceed by analyzing the station- 
ary amplitude-phase characteristics of the dot as a func- 
tion of the parameters v^c an d v a c, the dimensionless 
applied voltages. To proceed in this direction, we will 
assume that A and x vary slowly in time. Then, substi- 
tuting Eq. (6) into Eq. (2c) and the resulting expression 
in Eq. (2a), after averaging over the fast oscillations one 
obtains the following coupled differential equations for 
E(t) and x(r), 



dx 



X 



dE ■ 



(9a) 
(9b) 



Here, is the generating Hamiltonian function, 



- — I d01n{l + e WEcos\9)+v st +v ac cos{e- x )]y ( 9c ) 



FIG. 3: (Color online) Threshold ac voltage zy* c as a func- 
tion of the stationary dc voltage v st . Red bars correspond 
to results form numerical integration of Eqs. (2) while blue 
lines refer to the analytical perturbative scheme carried out 
in Eq. (10). The plot is calculated for a gold quantum dot of 
radius r = 4 nm and mass m = 5 x 10 -21 kg with uo — 10 GHz, 
Q = 1000 , d ~ 2 nm, A = 0.1 nm, T = 100 GHz, uj g = 10 GHz 
and T — 10 K. Consequently, e = 0.1, a = 837, Co — 1. 



Further analyses from Eq. (9) indicate that given a 
certain stationary dc voltage v st there exists a critical 
amplitude of the ac voltage above which the para- 
metric excitation of the mechanical vibration gives rise, 
see Fig. 3. This fact can be understood in the frame of 
a perturbative analysis for small E. In this regime, we 
find that the instability criterion due to the parametric 
resonance is given by 



J_ < sin(2x a t) 



J — 7T 



d0 cos(2<9) tanh( 
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(10) 

From Eq. (10) one can also notice that the right hand 
side of the inequality vanishes in the limit v st —> and 
the system can not be parametrically excited. In the 
instability criterion Xst is defined through the equation, 



djf(E, Xst ) 



dE 



= 0. 



(11) 
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By substituting Xst into Eq. (10) it is possible obtain 
a relation = ^aci^st) f° r the critical ac voltage as 
a function of the stationary dc voltage. This thresh- 
old ac voltage is shown in Fig. 3, where red bars corre- 
spond to results coming from the numerical integration 
of Eqs. (2) while blue lines refer to the analytical per- 
turbative scheme carried out in Eq. (10). The agreement 
between the two approaches is noticeable. 

In the remainder of this work, we will be mainly 
focusing on the case of exact resonance, Co = 1 + 
(€/4tt) J^d0[e- Ust + e^c cos(0- x )]-i 5 and would like to 

discuss in detail the outcomes from the stability anal- 
ysis. Due to the periodicity of the generating Hamil- 
tonian function, J^(E^x) — ^(^,X + ^J? the station- 
ary solutions of Eq. (9) come in pairs: to any solution 
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FIG. 4: (Color online) Stationary phase (a) and amplitude (b) 
of the system as a function of the applied dc voltage v s t — 2 
(orange square) and v s t — —2 (blue circles) for different dc 
voltages. The phase is almost zero for the positive dc voltage 
while it is non-vanishing for the negative dc voltage. Current 
transport is more feasible at negative dc voltages. The plot is 
calculated for e = 0.1, a = 837, Co — 1. 
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is displayed in Fig. 5. From the contour map, it is ev- 
ident that the conductive behavior of the system with 
respect to the stationary dc voltage v st is asymmetric. 
Therefore, in light of the result previously discussed, the 
values of the applied voltages should be choosen in or- 
der to maximize the charge transport. In this case, from 
Fig. 2, we look for a stationary phase %3 = ±tt/2 and, 
from Fig. 3(a), this condition is more likely achieved for 
negative values of v st , as shown in Fig. 5. 

To conclude, we have analized parametric effects in a 
nanoelectromechanical single electron device in the form 
of a three-terminal shuttle system under the influence of 
a gate voltage. In particular, we have shown that in the 
parametric excitation regime, i.e., when the frequency of 
the gate voltage is approximately equal to the mechanical 
frequency of the nano-oscillator, electrons are transferred 
between the electrodes maintained at the same potential. 
In this case, the direction of the shuttle current is con- 
trolled by the phase shift between the mechanical vibra- 
tions and gate voltage oscillations. 

All authors contributed equally to this work. Finan- 
cial support form the Swedish Research Council (VR) is 
gratefully acknowledged. 



FIG. 5: (Color online) Dimensionless current |/| obtained by 
numerical integration of Eqs.(2) as a function of the applied 
voltages. Current through the nanostructure is significal for 
negative ac voltage values. In the plot e = 0.1, a — 837, 
uj = 1. 

Pi — {XiiAi} corresponds a conjugated solution Pi = 
{Xi+ 7T , Ai}. Computational analyses have shown that for 
Vac(vst) < v* c (v s t) the system in Eq. (9) possesses four 
formal stationary points Pi = {7r/4, 0}, P^ = {37r/4,0} 
and their conjugates (for u st > 0, Pi (Pi) and P2CP2) are 
stable and unstable, respectively. For v st < 0, vicev- 
erse. If v a c{v s t) > v* c (v st ), two stationary stable points 
P3 = {X3? ^3} an d P3 appear on the phase diagram while 
the points Pi and P2 (and their conjugates) become un- 
stable. 

In Fig. 4, the phase shift and amplitude of the stable 
solution P3 = {x3,^3} are shown as a function of the 
applied ac voltage for two different dc voltages: v st = 2 
(blue circles) and v st = —2 (orange squares). From this 
graph one can notice that for the positive dc voltage, the 
phase shift and, as a consequence, the average current 
(see Eqs. (7) and (8)) are almost zero. Whereas for the 
negative dc voltage the phase is nearly tt/2 and the cur- 
rent increases approximately linearly with the amplitude, 
/ ^ As^euJo, where k is a numerical factor of order one. 

Finally, we investigated the current behavior as a func- 
tion of the applied voltages I = l(y s t-> v ac)- The result 
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